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Abstract 

The Casimir effect giving rise to an attractive force between the closely 
spaced two concentric spheres that confine the massless scalar field is cal- 
culated by using a direct mode summation with contour integration in 
the complex plane of eigenfrequencies. We devoleped a new approach ap- 
propriate for the calculation of the Casimir energy for spherical boundary 
conditions. The Casimir energy for a massless scalar field between the 
closely spaced two concentric spheres coincides with the Casimir energy 
of the parallel plates for a massless scalar field in the limit when the di- 
mensionless parameter 77, (77 = ^= where a (6) is inner (outer) radius of 
sphere), goes to zero. The efficiency of new approach is demonstrated by 
calculation of the Casimir energy for a massless scalar field between the 
closely spaced two concentric half spheres. 

PACS number(s): 03.70.+k, 12.20.DS, ll.10.Gh 



1 INTRODUCTION 



In 1948 Casimir p] had shown that two uncharged perfectly conducting 
parallel plates attracted each other. The main idea under this attractive force 
is that changes in infinite vacuum energy of the quantized electromagnetic field 
can be finite and observable [31 [3J 0]. The sign of Casimir energy depends 
on manifolds of different topology, geometry, the shapes and compositions of 
objects as well as the boundary and curvature [SJ [HI E] . The Casimir effect has 
been studied in different areas of theoretical physics and an important results 
may be found in detailed reviews [H [H [10] . 

It is well-known fact that the negative sign of the Casimir energy between 
the plates for the electromagnetic field produces an attractive force. Stimulated 
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by this event, Casimir hoped that one could construct an electron model as a 
perfectly conducting charged spherical shell. However, Boyer [TJJ showed that 
the quantum electromagnetic vacuum energy of a conducting spherical shell 
is positive which means that the force on the shell is outward and depends 
only on the radius of the spherical shell. This outward force was opposite 
to what Casimir had expected. Later, Milton et. al. [H] confirmed Boyer 
result by using Green's functions. More recently, the direct mode summation 
technique and the Casimir energy calculations for the spherical shell has been 
advanced by [131 [H]. Nowadays, the Casimir effect between two objects has 
also been considered [HI QH Q71 [H DH EUJ- Several different methods 
have been developed to calculate the Casimir effect beyond the proximity force 
approximation. One of the pioneering study is the work by M. Bordag [15] . 
where he developed to compute the small separation leading order terms of the 
Casimir interaction between a cylinder and a plane by using the path integral 
approach. Moreover, the exact result for the Casimir interactions between a 
finite number of compact objects of arbitrary shape and separation have been 
developed by using the functional determinant or multiple scattering approach 
[TBI H71 [T8l [T9l [201 [21] . Besides, the Casimir force problems of two concentric 
spheres [H [21 EH H3 I2FJ and tw0 concentric cylinders g3 [Ml HH] have also 
been calculated by using different regularization method recently. 

In this work we calculate the Casimir energy for the massless scalar field in 
a nontrivial two smoothly curved objects at close separation. The system we 
consider is made of two compact closely spaced spherical surfaces and the vac- 
uum gap consists of between two concentric spheres with radii a and b (b > a) . 
There are two motivations for us to perform the Casimir energy calculation for 
a massless scalar field between the closely spaced two concentric spheres, and 
between the closely spaced two concentric half spheres [2BJ [3D] ■ First it is in- 
teresting to find what similarities and difference of the Casimir energy of the 
toplogically similar geometry and comparing the parallel plates and the spher- 
ical shell, respectively. The sign of the Casimir energy will decide whether the 
Casimir force will be attractive or repulsive. We hope that there will be an 
application of our result to special systems in nanotechnologies and nanoelec- 
tromechanical devices. Second from the mathematical point of view the direct 
mode summation approach to the our geometry has been developed. Our aim is 
to show the simplicity and efficiency of the direct mode summation by contour 
integration when calculating the Casimir energy for a difficult boundary as the 
closely spaced two concentric spheres. The direct evaluation of the infinite sum 
over all the vacuum energy eigenvalues of the massless scalar field modes implies 
that the translation matrices in multiple scattering approach is not needed in 
order for the finite values for the vacuum energy to be obtained for given spher- 
ical boundary conditions. Moreover, the improved value for the computation of 
the Casimir energy for a massless scalar field between the closely spaced two 
concentric half spheres we reconsider here is developed by using the Abel-Plana 
sum formula for evenly spaced frequency spectrum for large argument. 

The organization of the paper is as follows. In Sec. 2. The Casimir energy of 
a massless scalar field subjected to spherical boundary conditions on between the 
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closely spaced two concentric spheres is calculated without any approaximation 
techniques. This approach will be employed for the Casimir energy between the 
closely spaced two concentric half spheres in Sec. 3. Concluding remarks and 
discussion of the Casimir energy for a massless scalar field in an annular region 
of our geometry is presented in section 4. 
The units are such that h = c = 1 . 



2 CASIMIR ENERGY BETWEEN TWO CON- 
CENTRIC SPHERES 

We start with the spacetime metric 

ds 2 = dt 2 - (dr 2 + r 2 d9 2 + r 2 sin 2 9d<t> 2 ) (1) 

in spherical coordinates. Where 9 <G [0,7r] and <j> G [0,27r]. The wave equation 
that the massless scalar field satisfies in these coordinates is given by 

□*(t,r,0,0)=O. (2) 

Where □ is the D'Alembertian operator associated with the metric given by the 
line element Eq. (1). Solution of equation (2) could be easily found by using 
the method of separation of variables and is given as 



(t, r, 0, <t>) = e ~ <W ' T ~ h [ A e Jt+i/2(ur) + Be N e+1/2 {ur)] Y £m (9, 4>) . 

(3) 

Where the Yi m (9, <f>) are the Spherical harmonics, and j£ + i/ 2 (wr) and N e+ i/ 2 (ujr) 
are the Bessel functions of the first and second kind. The coefficents Ai and 
Be can be determined from the boundary conditions, respectively. And £ = 
0,1,2,3,... and m = -£, -i + 1, ....0, 1, 2, ...,£ - l,£. 

We now impose the boundary conditions for our geometry (a (6)is inner 
(outer) radius of sphere) that is 



^utm (*, r = a,Q,4>) = 0, and similarly ^^im {t, r = b, 6, 4>) = . (4) 
The eigenfunction that satisfy the boundary conditions is 



* w , m (t, r, 9, 0) = c e- iw *r-5 



J 1+1/2 (wr) 



N e+1/2 {u>a] 



Ne +1/2 (ujr) 



Y tm (e,4>). 



(5) 

Where Co is the normalization constant and oj is the root of the following tran- 
scendental equation 



J e+ i/ 2 {ojb) N e+1/2 (uja) - Je +1/2 (uja) N e+1/2 (ujb) = 



(6) 
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The Casimir energy for the massless scalar field between two concentric 
spheres is 

^ oo £ oo 

Ec = 2 ^ 51 51 w ™ f ' 

= ^J^^J w n« where i/ = ^+-. (7) 

£=o n=i 

Where u;„£ are eigenfrequencies which are determined by solving the frequency 
equation given in Eq. (6). We need to describe the eigenfrequencies spectrum 
for the Casimir energy in a nontrivial smoothly curved objects at close sepa- 
ration. We note that the frequency equation involves an infinite series since 
the complete solution of Bessel's equation has series of ascending powers of lu. 
Bessel's series equation are convergent for all values of argument. But, when 
| a; | is large, the series converge slowly. Although the series has the converge 
slowly, the initial terms of such a series gives no information about the sum. 
Thus Bessel functions are needed to describe the transition for large argument. 
To overcome this difficulty, we need the rapidly convergent evaluation of the 
Bessel's function formula. The meaning of rapidly convergent is that the series 
rapidly approaches a constant, taking the limit as £ — > oo. To this aim, we 
will use the uniform asymptotic expansions of the Bessel's functions. The uni- 
form asymptotic expansions are useful in describing the transition of behaviour. 
Moreover, after taking the limit co —> oo at fixed £ the spectrum will consist 
of a discrete set embedded in a continum part. Hence, we should examine the 
behavior of the eigenfrequency spectrum for large arguments at fixed I. Thus, 
to carry out the summation with respect to £ in Ec, the sum ^n=i given 
in equation (7) replaced by X) n =l + S^Li <^ni where uj n i is the eigenvalue 
spectrum of the limit uj — > oo at fixed £ Then, Casimir energy which is 

defined by the eigenfrequency spectrum for large arguments at fixed £ and large 
order as £ — > oo can be written as 

oo oo oo oo 

E C = ^ V ^ U nl + ^ V ^ ^rd ■ (8) 

1=0 n=l 1=0 n=l 

Now, to calculate the eigenfrquencies for large arguments at fixed v, we use 
Hankel's asymptotic expansion |31) when v is fixed, Loa 3> 1 and cob 3> 1, we get 




(10) 
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and similar expressions for J v (cub) and N v (uib) with a interchanged for b. 
Putting (9) and (10) in the frequency equation given by (6), we obtain the 
zeros of frequency equation are almost evenly spaced 

u n i— I I H — ^ where n = 1, 2, 3, 4, 5, 6, ... . (11) 

The frequencies equation ui n g the first sum is given in Eq. (8) as the uniform 
asymptotic expansions of the Bessel functions at large £ as £ — > oo can be written 

as 

f v (vu, A) = J v (vw) N v (vw\) - J v (vbj\) N u (vu) (12) 

Where A = f (b > a). 

Then, the scalar Casimir energy between the closely spaced two concentric 
spheres can be written as 



oooo oooo/y \2 2 

t=0 n=l £=0 71=1 ' V ' 

= Ec + Ec . (13) 

Where u> n £ is the root of the frequency equation given in Eq. (12). 

We consider the first sum defined in Eq. (13). This divergent expression can 
be rendered finite by the use of a cutoff or convergence factor. Then we define 
the first sum, 



OO / 1 s oo 

= e( £+ £) (i4) 

where the factor of e~ au}nt plays the role of an exponential cutoff function, and 
Si = YlnLi Uni e~ au)nl is generated by the frequency equation (12). To evaluate 
the sum Si, we use the integral representation from the Cauchy's theorem [141 
that for two functions ft(z) and c/>(z) analytic within a closed contour C 
in which fe(z) has isolated zeros at x±.X2,X3 x n , 

±<f (fc^J^ln^)^^). (15) 
We choose <j>(z) — z e~ az where a is a real positive constant thus leads to 
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±-J cdz e- Z ^n Mz )=j: Zj e—. (16) 

Using this result to replace the sum Si by a contour integral, the first term of 
the Casimir energy becomes 

^ = §H) hl dze ~ aZz i l «^> (17) 

where the frequency function f v {vz) is given Eq. (12). The contour C encloses 
all the positive roots of the equation f v (vz) = 0. This contour can be conve- 
niently broken into three parts [T4j [26l [29] . These consist of a circular segment 
Cr and two straight line segments Ti and forming an angle <p an d tt — </> 
with respect to the imaginary axis. When the radius T is fixed, the contour Cr 
encloses a finite number of roots of the equation ] v {vz) = 0. Since the sum 
of these roots is obviously infinite, the radius T is a regularization parameter, 
taking the limit r —¥ 00 ( when a > ) means the removal of the regularization, 
the contribution of Cr vanishes provided that 0^0. Hence the exponential 
cutoff function in the Cauchy integral plays the role of the eliminate of the con- 
tribution to the circular part of the contour integral. Taking the contributions 
along r"i and T 2 which are complex conjugates of each other, then Eq. (16) 
becomes [H E] 



1 00 r°° d 

E c = lim V v Re e _i * / dy e - ioive y— In Uvye'^; a, b) . (18) 

n T^n Jo d V 



e=o 

Where 



Uiyye-*; a, b) = ~ [l v {vye^b) K v (vye^a) - I v {vye^a) K v {vye^b)] 

(19) 

Now we calculate the integral given in Eq. (18). Defining A = f (b > a) and 
after rescaling integral variable with yb —> y in equation Ec , one obtains 



1 f°° - d 

E c = r lim VfRe e _< * / dy e~ taye " P / b y— In Uvye'**, A) . (20) 

Where/^ye-^A) = -f [l v (vye^)K v (vye^\) - I^vye^XlK^vye-**)] 
stems from Eq. (19). 

Now we use the Lommel's expansions or the multiplication theorem for the 
function of f v {vye~^, A) [HI 132]. Thus, we have 
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fu{z,X) 



= I v (z) K v (zA) - I v (zA) K v (z) 



A -.y (A 2 ~l)' 



fe=0 



dz 



(21) 



Where A — 1 < 1. Applying the uniform with respect to z asymptotics for 
the modified Bessel functions at large v, after long calculations [25], we obtain 



vy — In f v (uye- l f,X) 
ay 



12 

(1-A 2 ) 4 



z/ye 



1-A 



2\3 



24 



2^5 



720 

1440 

(1-A 2 ) 6 
120960 

4 (vye-** 



(yye-**) 2 {-v 2 + 19) - {vye**)* 



3 (i/ye"*) 2 (^ 2 - 9) + 2 (vye-t+y 



[vye^) {Av A - 290v 2 + 1726) + (vye 



[Terms in even powers of {vye 



Inserting Eq. (22) into Eq. (20) and using the following integral result 



H2H) = e~^ I dye- iaye /b (uye~^) 

2n+l 



2)1 





n+1 



= i {-y^ 1 (2n)! v 
then Eq. (20) becomes 



7 \ * 

■In I M 

a) 



where n = 0, 1,2,3... (23) 
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E c = — 7 limV Re 

( ' 

fl-A 2 



oo 



£=0 
i2 n4 r 



a 



6 



+ i 



a 



12 



+i 



+i 



+i 



720 

:i-A 2 r 

720 

'1-A 2 ) £ 



2(-^ 4 + 19i/) - 



-6(^ 4 - 9^ 2 



24i/* [ - 



48^ 4 



580^ 4 + 3452zy 2 ) 



24(8j/ 6 - 149i/ 4 ) 



120960 

[Terms in imaginary number and even powers of v]} 



All terms in the above equation have the singular term in the regulator param- 
eter a and purely imaginary. Taking the real part of the parenthesis, thus it 
leaves the zero result. 



E c = 



(25) 



The meaning of this result is that there is no contribution from £ — ¥ oo modes 
for the Casimir energy between two concentric spheres [26 . Now we return to 
the second sum given in Eq. (13), included high eigenfrequency modes i.e. 
uj — > oo at fixed I. 



OO OO 

e=o ri=i 



b-. 



ab 



(26) 



Where v = (£ + |) . This divergent sum can be regularized by using the Abel- 
Plana sum formula. Before we proceed, we recall the Abel-Plana formula gives 
an expression for the difference between the sum and corresponding integral, 
which could be given as [9l [TO] 



2880z/ 6 



a 



(24) 



Reg 



.71=1 



71 = 



£j»- f(x)dx = -/(0) 



' ' e 2.t _ i ( 27 ) 



Where / (z) is an analytic function in the right half plane and Reg refers to the 
regularized value of the sum. The other useful Abel-Plana sum formula for the 
half integer number is 
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CO oo oo 

n-0 Q Q 

We can rewrite the sum given in Eq. (26) using the Abel-Plana sum formula, 
which leads to 



Er 



where £ 



= ^1/ Reg 



e=o 



E 



6- 



OO 

^ OO oo „ 

^E - 2 - 2 E- / 

^ ab i=o t=o i 



ab 



tir 



ir- 
ah 



— 



dt 



o2TTt 



- 1 



(29) 



2d 



and d = b — a. The first divergent sum in Eq. (29) can be 



removed by using the Hurwitz zeta function [33] 



E 

£=0 



Thus Eq. (29) becomes 



(30) 



Er 



where F (v) 



= -2£ F(y) 



1=0 



tn 

b — a 



2 2 



ab 



dt 



27Tt 



(31) 



Using the half integer Abel-Plana sum formula given in Eq. (28), one obtains 



Er 



1 



47T d 3 



-C(4) + 2 ? : 



F(it) - F(-it) 



dt. 



(32) 



Where £ (s) is the Riemann zeta function. The divergent term occurs in the 
second term of the Eq. (32). We avoid this problem by substituting v — i € for 
v with the understanding that we will let €— > at the end. Thus the Casimir 
energy for a massless scalar field between the closely spaced two concentric 
spheres is obtained 
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E, 



c 



E ( 



Ec 



1 



47T d 3 



-C(4) 



.C(2) 



1 



12' C(4) 



= 7r " q6 
360 d 3 



I , JLiL 
47r 2 ab 



(33) 



Where 77 = ^7=. And known formula £ (4) = §5 and ((2) = ^1. The main 
contribution to the Casimir energy is given by the second term which represents 
the high frequency modes at fixed I between the closely spaced two concentric 
spheres for a massless scalar field. The Casimir energy per unit surface area on 
the inner sphere (the total surface area A = 47ra 2 ) can be written as 



Ec 
A 



16tt 2 



C(4) 
d 3 



1 H 77' 

12 ' 



:C(2) 

C(4) 



(34) 



This result is interest at the limiting case which is narrow slit is defined by 



V 

\/~ab 



y/ab 



< 1 



We easily analysis that in the limit b — > a (i.e. 7/ — > and 



^2 — > 1) which means that the close separation between two spheres one finds 
that the leading term of the Casimir energy per unit area can be written as 



Ec 
A 



1 C(4) 



(35) 



16tt 2 d 3 ' 

This result is exactly the same as the Casimir energy of the parallel plates for 
a massless scalar field [33]. Thus our approach developed here has been the 
satisfactory check. 



3 CASIMIR ENERGY BETWEEN TWO CON- 
CENTRIC HALF SPHERES 

The our approach can be easily employed to the computation of the Casimir 
energy of a massless scalar field between the closely spaced two concentric half 
spheres. In this case of our spherical boundary }26j the sum in Eq. (13) can be 
written as 

*-iE«E-*+5£'£V(£:) + 3- (36) 

l=\ n=l l=\ n=l » V / 

Where v = £ + i and uj n g is the root of the transcendental equation given in 
Eq. (6). From Eq. (15) we obtain 
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E c — E c + E c 

1 °° / 1 \ Z" 00 --<4 d 

£=1 v 7 n=l V \ / 

Where In f v (yye~ % ^, A) given in cq. (22). The first term (= i?c) is equal 
to zero from the Eq. (25). Thus the Casimir energy between the closely spaced 
two concentric half spheres becomes 

Using the Abel-plana sum formula one obtains 



1 OO j oo 

E C = -J-Y(S-%)-±YF(v). 



(39) 



OO 

where F (u) = ~ y) / (j/ 2 -!)' 



/2 



e 27T !/ _ i ' 



and d = b — a. 



Again using the half integer Abel-Plana sum formula for the second sum in Eq. 
(39), we obtain 



E c 



1 



8tt <P 



C(4) 



1 T] 



7T C(3) 7T 



1 



4''C(4) ' 24 V C(4) ' 12'' C(4) 4tt'' £(4) 

(40) 

We have used the Hurwitz zeta function i.e. Q{—m,q) = — ^^fi^ ' m = 
0, 1, 2, 3, ... where B m+ \ (q) is the Bernoulli polynomials. 

The Casimir energy of a massless scalar field between the closely spaced two 
concentric half spheres per unit half surface area on the inner sphere (the half 
surface area A = 2ira 2 ) can be written as 



1 

777*7 



2 C(2) 



1 3 1 

— rf 



Ec 



1 

'167T 



C(4) 

d 3 



^(3) 7T o 1 



+ 777?? 



+ 777?? 



1 „2<(2) 



1 3 1 



4'C(4) 24' C(4) 12' C(4) 4tt ' C(4). 

(41) 

Where £ (s) is the Riemann zeta function. 
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Taking the small separation limit (b — > a and 77 — > 0) one obtains 

E c 1 C(4) 



(42) 



A 16tt d 3 ' 

This result coincidences with the Casimir energy of the parallel plates for the 
massless scalar field [Ml. 



4 CONCLUSION 



In the present paper we considered the quantum vacuum energy for a mass- 
less scalar field between the closely spaced two concentric spheres, and between 
two concentric half spheres at small separations. The annular region in our 
geometry was considered since all massless scalar fields in the region r < a and 
r > b are equal to zero. The Casimir energy for a massless scalar field in annu- 
lar region was evaluated by a direct mode summation method. We have used 
the explicit expression for the frequency equations which include the product of 
Bessel functions. In order to evaluate the product of Bessel function expansion 
for large order in our approach, the Lommel's expansions has been used in the 
Casimir problem for the first time. These calculation is a direct application of 
the principle of argument from the complex integral with the cutoff function, in 
direct analogy to the Casimir calculation, together with a modification of the 
contour in order to ensure the convergence of the Cauchy integral expressions. 
Moreover, the product of Bessel function expansion for the large argument is 
defined by evenly spaced eigenfrequency spectrum. Although we use the Abel- 
Plana sum formula, divergent terms appear in our calculations. To remove this 
divergence we have applied the formal technique of the Zeta function regular- 
ization. Thus, our calculation implies that further regularization is not needed 
in order for the finite values for the vacuum energy to be obtained for given 
boundary conditions. 

We are particularly interested in calculating the Casimir energy between 
spherical surfaces close to each other is that any approximation technique is not 
needed. The interesting point of our calculations is that all contributions in 
the Casimir energy for a massless scalar field comes from the higher frequencies 
for fixed I between two surfaces boundary conditions. £ — > 00 frequency modes 
contribution in the Casimir energy is zero for close separation of annular region 
in present geometries. 

Although the Casimir energy sign for the massless scalar field on spherical 
shell with the Dirichlet boundary condition is positive, our analysis reveals the 
negative sign of the Casimir energy between the closely spaced two concentric 
spheres, and between two concentric half spheres at small separations. Then, 
the Casimir force for a massless scalar field between two concentric spheres, and 
between two concentric half spheres are always attractive which is the same per 
unit area on a pair of parallel plates for a massless scalar field [31] . Moreover, the 
leading term of the Casimir energy between the closely spaced two concentric 
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spheres given in Eq. (33) agrees with that obtained the Casimir effect between 
two spheres at small separation by using the functional determinant or multiple 
scattering approach [2TJ1 I3"5] . 

Closing, it is worth noting that, as far as we know, such boundary conditions 
with between the closely spaced two concentric spheres has been considered in 
the Casimir problem with the use mode sum technique for the first time. As far 
as we know this result is obtained here for the first time. 

For future works, it would be interesting to consider the electrodynamics 
Casimir energy between the closely spaced two concentric spheres with those 
obtained. 
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